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RESUMEN 

We discuss the role of the electron temperature in abundance determinations in ionized nebulae (planetary 
nebulae and giant H n regions). We show that, even when observations provide a direct estimate of T e , abun- 
dance determinations may sometimes be significantly in error for reasons other than hypothetical temperature 
fluctuations or uncertainties in ionization correction factors, and we show striking examples both in the high 
and the low metallicity regimes. On the other hand, even without a direct measure of T c and under certain 
conditions, it is possible to give rough estimates (or limits) on the abundances, both in GHRs and in PNe. 
Although the main mechanisms determining the electron temperature in ionized nebulae have been known for 
half a century, we emphasize that there are still important problems to be solved. Finally, we advocate for a 
different description of temperature inhomogeneities than the scheme generally used. 

ABSTRACT 

We discuss the role of the electron temperature in abundance determinations in ionized nebulae (planetary 
nebulae and giant H n regions). We show that, even when observations provide a direct estimate of T e , abun- 
dance determinations may sometimes be significantly in error for reasons other than hypothetical temperature 
fluctuations or uncertainties in ionization correction factors, and we show striking examples both in the high 
and the low metallicity regimes. On the other hand, even without a direct measure of T c and under certain 
conditions, it is possible to give rough estimates (or limits) on the abundances, both in GHRs and in PNe. 
Although the main mechanisms determining the electron temperaturein ionized nebulae have been known for 
half a century, we emphasize that there are still important problems to be solved. Finally, we advocate for a 
different description of temperature inhomogeneities than the scheme generally used. 
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1. INTRODUCTION 

Emission lines, and especially the forbidden lines 
used for abundance determinations in ionized neb- 
ulae are strongly sensitive to the electron tempera- 
ture, T e . It is common thought that, when T e can 
be derived directly from observed spectra, for exam- 
ple using the [O m] A4363/5007 ratio, abundance 
determinations are fairly reliable. We will show that 
this is not necessarily the case (§3). We will also 
discuss the reliability of abundance determinations 
when T e is not measured directly, both in the case 
of planetary nebulae (PNe) and giant H n regions 
(GHRs) (§4). Then, we will address the question of 
whether temperatures observed in ionized nebulae 
are compatible with the predictions from photoion- 
ization models and we will touch upon the problem 
of temperature inhomogeneities (§5). 

2. BASICS 

As is known from general considerations on the 
thermal balance in ionized nebulae (Spitzer 1948, 



1949, Spitzer k Savedoff 1950, Osterbrock 1989), the 
electron tempcratureresults from a balance between 
energy gains and losses. It is easy to show that, 
schematically, the energy gains per unit volume per 
unit time in a photoionized nebula can be written 
as: 

G = n 2 e a A (H°,T e )^E 

where n e is the electron density, otA(H°, T) is the 
recombination coefficient of hydrogen, and E is the 
mean energy of the absorbed photoelectrons, which 
is roughly rooportional to the effective temperature 
of the ionizing radiation field, T*. Energy losses are 
due to a variety of processes, emission of collisionally 
excited lines being the most important. In a two- 
level approximation and neglecting collisional deex- 
citation, they can be schematically written as: 

L = Yn e n(X«)^T-°- 5 e-^ kT °his l 
where n(A y ) is the number density of ion X l \ and 
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the index / refers to a specific line, f2; being the col- 
lision strenght, and Xl the excitation level. 

Thus, the electron temperature in a nebula is not 
a strong function of the distance to the ionizing stars, 
since it does not depend on the flux of ionizing pho- 
tons. It is larger for higher values of T* and lower 
values of the metallicity, Z. Fig. 1 illustrates the 
heating and cooling rates in a photoionizcd nebula, 
at two different densities. At high density, the in- 
frared fine-structure lines are partly quenched by col- 
lisional deexcitation, resulting in a larger T . Note 
also that, because of the shape of the cooling curves, 
there is a large temperature difference between the 
high and the low excitation zones at high metallici- 
ties. 




Fig. 1. Approximate heating and cooling rates divided by 
nl , in arbitrary units. Circles: approximate cooling rates 
in the high excitation zone (0 ++ ions only). Triangles: 
cooling rates in the low excitation zones (0 + and N + 
ions only). Left: n e = 100 cm -3 ; right: n e = 10000 cm -3 . 
The heating rate is represented by crosses. By shifting 
the heating curve with respect to the cooling curves, one 
can understand how T c varies with metallicity or with 
the effective temperature of the ionizing radiation field. 



3. WHEN T e -BASED ABUNDANCE 
DETERMINATIONS ARE IN ERROR 
3.1. One example at low metallicity: the effect of 
collisional excitation of Ha 
At low metallicities, because of the high T , col- 
lisional excitation of the Ha line becomes non negli- 
gible and the Ha/H/3 ratio no longer has the recom- 
bination value. If the dereddening procedure of an 
observed spectrum involves Ha (which is commonly 
the case), this translates into an overestimation of 
the extinction C(H/3), and an error in the dereddened 
line ratios. It can easily be shown that this error in 
the line ratios is independent of extinction. It can be 
large when the lines in question have very different 
wavelenghts. 
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Fig. 2. Photoionization models for evolving PNe with Z 
= 0.1 Zg, M ne 6=0.2 Mq, v exp =20 km s~ . a) computed 
Ha/H/3; b) case B Ha/H/3; c) C ill] A1909/[O in] A5007 
"dereddened" using case B with respect to actual C ill] 
A1909/[O in] A5007; d) same as c) for He I A5876/H^. 



This is illustrated in Fig. 2, which shows the 
results of photoionization models at metallicity Z— 
0.1 Zq for evolutionary sequences of planetary neb- 
ulae whose central stars follow the theoretical tracks 
of Blocker (1995). These models were computed us- 
ing the procedure outlined in Stasihska et al (1998) 
Q. We see that, at such a metallicity, there is quite 
a range in epochs during the evolution of PNe (ma- 
terialized in the figure by the surface brightness in 
H/3), where Ha/H/3 is significantly larger than the 
recombination value. The resulting error in the de- 
rived helium abundance may reach 5 - 10%. A sim- 
ilar effect is expected in very low metallicity H n 
regions. This is important to keep in mind for the 
determination of the pregalactic helium abundance 
from low metallicity GHRs. Davidson & Kinman 
(1985) indeed warned about this effect, but it does 
not seem to have been properly taken into account 
in later studies (see e.g. Olive et al. 1997 or Izo- 
tov et al. 1999 and references therein). The error 
in C/O obtained using the C in] A1909 and [O III] 
A5007 lines may reach almost a factor 2 in the PN 
models shown here. Again, a similar effect is ex- 
pected in low metallicity GHRs. This could explain 



1 Unless otherwise stated, the photoionization calcula- 
tions presented here used the code PHOTO as described by 
Stasihska & Leitherer (1996). 
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why the carbon abundance estimated in I Zw 18 is 
larger than expected for its metallicity (Garnett et al 
1997a). However, in the case of UV data, a possibly 
even more important problem is the uncertainty in 
the reddening law (Garnett et al 1999). 

3.2. One example at high metallicity: the absence of 
[O ill] optical lines from the ++ zone 

At high Z, because of the intense cooling by 
the forbidden lines, T c is low (see Fig. 1) and 
the transauroral lines used for T c diagnostics be- 
come hard to detect. So far, there has been no 
T e -based measurement of metallicities higher than 
about 0.5 Z Q in any GHR. With 8 m-class tele- 
scopes, however, it becomes possible to detect weak 
[O in] A4363 fluxes, and the hope has been ex- 
pressed that it will be possible to directly probe 
oxygen abundances at metallicities higher than so- 
lar. However, as already shown by Stasihska (1978) 
and Garnett (1992), high metallicities induce an im- 
portant temperature drop in the ++ zone because 
cooling is dominated by the [OIII]52^ and [OIII]88/x 
lines which are almost independent of T (see Fig. 
1). Therefore, the value of T derived from [O in] 
A4363/5007 would largely overestimate the temper- 
ature characteristic of the ++ zone. Consequently, 
classical T e -based methods would strongly underes- 
timate the oxygen abundance. This is illustrated in 
Fig. 3. The derived O/H ratio always stays below 
the solar value, and would be underestimated by or- 
ders of magnitudes for metallicities over twice solar 
if such methods were applied. Incidentally, note the 
rise in the [O in] A4363 flux at log O/H + 12 > 9.1. 
This is due to the fact that recombination becomes 
the dominant process for the excitation of this line 
at the low temperatures involved. 

Does this mean that it is impossible to reliably 
determine the oxygen abundance in the central parts 
of galaxies, where it is supposed to be oversolar? The 
most straightforward way, of course, would be to use 
the [OIII]52/x line (Simpson et al. 1995, Rudolph 
et al. 1997). Another possibility, suggested by 
Stasihska (1980), would be through high spatial reso- 
lution emission line imaging. Metal rich GHRs would 
show very characteristic and unusual spatial correla- 
tions between line intensities, since [O in] and [O ii] 
and would be emitted almost cospatially in the high 
T c zones, while the hydrogen and helium recombi- 
nation lines would preferentially arise in the low T e 
zones. Of course, detailed photoionization modelling 
reproducing the observed tendencies would still be 
necessary to provide a quantitative estimate of O/H. 

Planetary nebulae provide an interesting substi- 
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Fig. 3. Results from a simple model for a GHR having 
T*=45000 K, Q H =3 10 51 ph s"\ n=30 cm" 3 , filling 
factor = 0.1. Left: predicted flux in the [O ill] A4363 
line for an object at a distance of 1 Mpc; right: O/H 
derived by the classical T e -based method as a function of 
input O/H. 
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Fig. 4. Same as Fig. 3 for a simple model ofa ionization 
bounded PN having T*=150 000 K, Q H =3 10 47 ph s _1 , 
n=10 4 cm~ 3 , filling factor = 0.3 



tute to GHRs for abundance determinations in the 
central parts of galaxies. As shown by Richer et al. 
(1998), the most luminous PNe are good tracers of 
the ISM oxygen abundance. Such PNe, being young 
(therefore dense) and excited by stars with T* of 
the order of 100000 K, have larger T c than GHRs of 
same metallicity (see Fig. 1). Therefore, they can be 
detected in [O in] A4363, even though their intrin- 
sic luminosity is orders of magnitudes smaller than 
that of GHRs (compare Figs. 4 and 3). In addition, 
electron temperature gradients are not expected to 
be so severe as in GHRs (as can be seen from Fig. 
1), so that the classical T c -based method to derive 
O/H should give reasonable results. Indeed, Fig. 4 
shows that, at least theoretically, this is the case. 

4. ESTIMATING O/H WITHOUT T 
4.1. Strong line methods for GHRs and PNe 
It has been shown by Pagel et al. (1979) that, 
in GHRs, provided that the hardness of the ionizing 
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radiation field is related to the metallicity, it should 
be possible to calibrate some well chosen strong line 
ratios as a function of O/H (see Stasihska 2001 for 
a summary of the basics and a short bibliography of 
strong line methods). The extremely tight correla- 
tion between ([O m] A5007 + [O n] A3727) / H/3 and 
T -bascd O/H determinations in high signal-to- noise 
observations of H II galaxies (see Fig. 5) is an em- 
pirical confirmation that this method works at low 
metallicities, as was suggested by Skillman (1998). 
Incidentally, this implies that GHRs in blue com- 
pact galaxies constitute a very homogeneous class as 
regards their excitation and nebular properties (see 
discussion in Stasihska ct al 2001). 
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Fig. 5. ([O in] A5007 + [O n] A3727) / H/3 versus O/H 
for GHRs in a sample of blue compact galaxies with T c - 
based derivations of O/H. Data from Izotov & Thuan 
(1999 and references therein). The number of objects is 
indicated in parenthesis. 
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PNc? Since PN central stars span a wide range 
in T* (30 000 K - 200000 K), one does not expect 
such methods to work in the same way as for GHRs. 
However, Fig. 6 shows that PNe outline a rather 
well defined upper envelope in the ([O in] A5007 + 
[O II] A3727) / H/3 vs. O/H plane, as already noted 
by Richer (1993) for Magellanic Clouds PNe. This 
upper envelope probably corresponds to the hottest 
central stars, and can be used to infer a lower limit 
for O/H. Such a procedure has been used by Richer 
& Mc Call (1995) and Richer et al. (1998) for extra- 
galactic PNe with [O in] A4363 undetected. 

4.2. Photoionization modelling 

Contrary to a widespread opinion, photoioniza- 
tion models (be it of PNe or of GHRs) are usually 
incapable of providing abundances if T e is not di- 
rectly measured. It is easy to produce models dif- 
fering only in abundance ratios that will reproduce 
all the observed emission lines. This is illustrated 
in Fig. 7, which shows the run of [O in] /H/3 and 
[O n]/H/3 as a function of O/H for a simple PN 
model. The reason of such a behavior is the same for 
PNe and GHRs and has been wideley commented in 
the literature. As O/H increases, T c decreases due 
to enhanced cooling, and the optical forbidden lines 
are gradually overtaken by the infrared lines, lead- 
ing to these bell-shaped curves. The maximum of 
the curves occurs at a metallicity around solar for 
dense PNe ionized by stars with T* > 100 000 K, 
and around half solar for GHRs. 
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Fig. 7. Variation of [O in] /H/3 and [O n]/H/3 as a func- 
tion of O/H for a simple density bounded model of PN 
with M ne6 =0.5 Mq, r*=100 000 K, Q H =3 10 47 ph s _1 , 
n=10 3 cm' 3 . 
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Fig. 6. Same as Fig. 5 for PNe in various galaxies: M31, 
the bulge of the Milky Way, the LMC and the SMC. Data 
as collected in Stasihska et al. (1998). 

Can such a strong line method be applied to 



Table 1 shows the results of photoionization mod- 
elling of the Galactic Bulge PN PNG 351.2+05.2 
(M 2-5) for which Ratag (1991, 1997) had claimed 
an abundance of one fourth solar. Using the Lviv 
photoionization code, Stasihska, Malkov & Golovaty 
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(1995, unpublished) produced two classes of models 
for this object accounting satisfactorily for the ob- 
served data (line strenghts, total flux and radius). 
The first one indeed has an oxygen abundance sim- 
ilar to Ratag's, but the second has O/H larger than 
solar! 

TABLE 1 
MODELS FOR THE PN M 2-5 



obs Ratag SMG SMG 

low Z high Z 

T* (K) 37500 39000 39000 

r* (10 10 cm) 5.0 4.9 

R in (pc) 0.062 0.065 

R OM t (pc) 0.10 0.087 0.085 

-logF(H/3) 11.409 11.407 11.411 

n (cm" 3 ) 2050 1800 1800 

He 0.117 0.117 0.100 



c 






1.50(-3 


1.20 (-3 


N 




4.80(-4 


2.50(-4 


6. 00 (-4 


O 




2.20(-4 


2.40(-4 


1.20(-3 


Ne 






5.00(-5 


2.40(-4 


S 




2.30(-5 


3.00(-6 


7.00(-6 


Oil 3727 


0.596 




0.587 


0.604 


Ne 3869 






0.014 


0.0096 


OIII 4363 


<.0013 




0.0006 


0.0001 


Hell 4686 






0.0004 


0.0003 


HI 4861 


1.00 




1.00 


1.00 


OIII 5007 


0.283 




0.304 


0.275 


NI 5200 


0.0149 




0.0043 


0.0087 


Nil 5755 


0.0071: 




0.0151 


0.0060 


Hcl 5876 


0.128 




0.126 


0.128 


OI 6300 






0.0054 


0.0116 


Nil 6584 


2.85 




2.79 


2.81 


SII 6717 


0.0565 




0.053 


0.0558 


SII 6731 


0.084 




0.077 


0.0826 


Oil 7325 


0.0091: 




0.0126 


0.0063 



This demonstrates that, in absence of a direct 
measure of T e , photoionization modelling by itself 
does not necessarily give a clear clue to the metallic- 
ity of nebulae when only optical data are available. 
In some cases, however, when solutions on either the 
low or high metallicity side can be rejected using 
arguments of astrophysical nature, photoionization 
modelling can provide estimates of O/H. 



An example where photoionization modelling 
gives interesting limites on O/H is the recently dis- 
covered PN in the galactic halo, SBS 1150+599A 
(now PN G135.9+55.9, for full details sec Tovmas- 
sian et al., 2001 and these proceedings). In that 
case, solutions at high metallicity can obviously be 
discarded. There is still a wide range of solutions 
compatible with the observed line ratios, but when 
additional considerations are taken into account, it is 
found that the object has O/H of the order 1/500 so- 
lar (making it possibly the most oxygen-poor object 
known in the Galaxy). 

5. DO WE UNDERSTAND THE 
TEMPERATURE IN IONIZED NEBULAE? 

So far, we have not questionned the electron tem- 
perature itself. As a matter of fact, there are indi- 
cations that simple photoionization models are inca- 
pable of reproducing satisfactorily the observed T e 
in quite a variety of situations. There are cases of 
GHRs as well as of PNe, where photoionization mod- 
els give too low a T c (Campbell 1990, Garcia- Vargas 
et al. 1997, Stasihska & Schaercr 1999, Luridiana 
et al. 1999, Luridiana & Peimbert 2001, Peha et al 
1998), calling for additional heating sources. In other 
cases, however, no T e problem is apparent (Gonzalez 
Delgado & Perez 2000) and in some cases, models 
seem to give too high a T e (Oey et al. 2000). So far, 
we are still lacking a synthetic view of these anoma- 
lies. One difficulty is that, even if a model reproduces 
the observed line ratios satisfactorily, this does not 
necessarily mean that the model is a fair representa- 
tion of the object. For example, the first attempts to 
solve the T e problem in I Zw 18 were not supported 
by later observations (see discussion in Stasihska & 
Schaerer 1999). 

It is of prime importance to understand what de- 
termines T e in real nebulae if one wishes to derive 
reliable abundances, since one way or another, the 
observed line ratios used to derive abundances from 
optical or ultraviolet data depend on T e . 

In the following, we comment on two disputed 
topics. 

5.1. Temperature gradients in photoionized nebulae 
While the electron temperature, being deter- 
mined by the hardness of the ionizing photons and 
not on their flux (c.f. §1), does not necessarily fall 
off with radius, some T c gradients are nevertheless 
expected. A consideration of the heating and cool- 
ing curves (Fig. 1) shows that T c should rise out- 
wards if the outer zones contain ions that are less 
efficient for the cooling. As mentioned above, the 
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effect can be quite strong at metallicity above solar. 
At low metallicity, on the other hand, T c decreases 
outwards. The hardening of the absorbed photons 
towards the ionization front, due to the frequency 
dependence of the photoabsorption cross-section is 
another factor that affects the radial T e distribution. 
The predicted effect is however rather mild (see e.g. 
the GHRs models of Stasinska & Leitherer 1996, or 
the grid of PN models from Stasinska et al. 1998). 

Temperature gradients can be estimated from the 
comparison of T c derived from various line ratios, 
such as [O in] A4363/5007and [N n] A5755/6584 (see 
eg. Garnett et al. 1997b for GHRs and Mc Kenna 
et al. 1996, Kinsburgh & Barlow 1994 for PNe). 
Recently, Pena et al. (2001) found that in PNe 
with Wolf-Rayet central stars, T e ([0 lll])/T e ([N n]) 
could be significantly different from 1 (by as much as 
40%) and that it seemed to decrease with increasing 
++ /0 + . Such a tendency cannot be reproduced by 
classical photoionization models. Dopita & Suther- 
land (2000) have shown that, by including the pho- 
toelectric effect on dust grains, one obtains larger 
T e gradients than in dust free models, because dust 
heating affects essentially the inner zones. The cases 
with T c ([0 iii])/T c ([N II]) < 1 observed by Peha et 
al. (2001) however remain unexplained. One possi- 
bility is to advocate a contribution from shock heat- 
ing in the outer zones, but quantitative models have 
still to be worked out. 



emitting regions, the mean electron temperature de- 
fined by Pcimbert (1967) becomes 



T n = 



Ti + fT 2 



1 + / 

and the fluctuation parameter becomes 



t 2 = (Ti~ To) 2 + f(T 2 - T ) 2 
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Fig. 8. Results from a two zone toy model defined by 
T =8 000 K and t 2 = 0.04 (see text). 



5.2. Temperature inhomogeneities 

In a seminal paper, temperature inhomogeneities 
have been advocated by Peimbert (1967) to explain 
the discrepancies between T as determined by var- 
ious methods. Peimbert introduced the tempera- 
ture fluctuation scheme and observational work ac- 
cumulated over the years by Silvia Torres-Peimbert, 
Manuel Peimbert and their coworkers and disciples 
seems to indicate a value of t 2 ~ 0.04. There are 
many questions about the reality, the ubiquity and 
the causes of these T e fluctuations (see e.g. papers 
by Esteban, Pcimbert, Binctte, Liu and Pequignot 
in this volume, and references therein). 

In fact, the temperature fluctuation scheme may 
not be quite appropriate to describe temperature in- 
homogeneities in ionized nebulae. A different ap- 
proach was presented by Mathis et al. (1998). Here, 
and in order to help visualising the problem, we use 
a toy model consisting of two homogeneous zones of 
volumes V\ and V2 with temperatures T\ and T2, 
electron densities n\ and 112, and densities of the 
emitting ions (e.g. ++ ) Ni and N 2 . If we call / 
the ratio (A^^V^/X-WiniVi) of the weights of the 
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Fig. 9. Results from a two zone toy model denned by 
T =15 000 K and t 2 = 0.04 (see text). 



In Figs. 8 (resp. 9) we explore the values of 
/ that lead to the canonical value t 2 = 0.04 at 
T o =8 000 K (resp. T =15 000 K), and show the 
effect on the ++ abundances calculated by using 
the T c derived from [O ill] A4363/5007. The case 
/ = 1, i.e. regions of equal weight, corresponds to 
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Ti=12 000 K and T 2 =8 000 K. This is quite a large 
temperature difference, which would correspond to 
heating- or cooling rates differing by about a factor 
3 between the two zones (see Fig. 1). No wonder 
that photoionization models have difficulties in pro- 
ducing such values of t 2 around these temperatures! 
When / ^> 1, there is a high weight zone at T 2 < T 
and a low weight zone at T\ ^ Tq. Such a situa- 
tion could correspond to a photoionized nebula with 
small regions being heted by shocks or conduction. 
When / <C 1, there is a high weight zone at T\ > T 
and a low weight zone at T 2 <C T which could cor- 
respond to high metallicity clumps. 

As expected, Figs. 8 and 9 show that the ++ 
derived from infrared fine structure lines is correct. 
The optical recombination line O II A4651 also give 
a correct ++ . However, high temperature dielec- 
tronic recombinations are not included in the present 
day atomic calculations. They might well boost O II 
A4651 in zones of high T c (around 20 000-50000 K, 
seeLiuetal. 2000, Dinerstein et al. 2000). This is an 
attractive explanation to the fact that in PNe, ++ 
abundances derived from recombination lines are far 
greater than derived from optical or infrared lines 
(see Liu, these proceedings). Compared to the al- 
ternate explanation invoking oxygen-rich condensa- 
tions (Liu et al 2000, Pequignot, these proceedings), 
this one has the advantage of being naturally ful- 
filled if additional heating mechanisms substantially 
raise the electron temperature in limited volumes of 
the emitting gas, such as shock or conduction fronts. 
Besides, they do not require to find a mechanism to 
produce oxygen-rich pockets in planetary nebulae. 

Following expectations, ++ derived from [O in] 
A5007 is generally underestimated, but it is interest- 
ing to note that the magnitude of the effect depends 
both / and on To. It is negligible when log / > 3-4, 
because [O in] A4363 saturates above ~ 50 000K. In 
the To = 15 000 K case, the effect is never strong 
anyway. For the To = 8 000K case, errors in ++ 
of up to a factor 2-3 can occur in the regime where 
[O III] A4363is significantly emitted in both zones. 
As for O ill] A1661, it may lead either to an underes- 
timation or to an overestimation of ++ /H, because 
the dynamical range of the line is larger than that of 
[O in] A4363. 

Fig. 8 and 9 thus amply demonstrate that the 
classical temperature fluctuation scheme can be mis- 
leading. Even in a simple two zone model, the situ- 
ation needs 3 parameteres to be described, not 2. In 
our representation, these parameters would be T\, 
T 2 and /, but other definitions can be used. This 
means that care must be taken when drawing con- 



clusions on temperature inhomogencities by compar- 
ing abundances derived from various lines. We argue 
that such diagrams as shown here might be useful to 
pinpoint the physical cause of temperature inhomo- 
geneities in ionized nebulae. 

This wandering - and wondering - about the elec- 
tron temperature in ionized nebulae is dedicated to 
Silvia Torres-Peimbert and Manuel Peimbert, at the 
occasion of a meeting in their honor. I am grate- 
ful to the Organizing Committee of for having made 
possible my participation to this stirry conference. 

REFERENCES 

Blocker, T., 1995, A&A, 299, 755 

Campbell, A., 1990, ApJ, 362, 100 

Davidson, K., Kinman, T. D., 1985, ApJS, 58, 321 

Dinerstein, H.L., Lafon,C.E., Garnett,D.R., 2000, in 
"Asymmetrical Planetary Nebulae II: From Origins 
to Microstructures" , ASP Conference Series, Vol. 199. 
Edited by J. H. Kastner, N. Soker, and S. Rappaport 

Dopita, M. A., Sutherland, R. S., 200, ApJ, 539, 742 

Garcia- Vargas, M. L., Gonzalez Delgado, R., Perez, E., 
Alloin, D., Terlevich, E., 1997, ApJ, 478, 112 

Garnett, D. R., 1992, AJ, 103, 1330 

Garnett, D. R., Shields, G. A., Skillman, E. D., Sagan, 
S. P., Dufour, R. J., 1997b, ApJ, 489, 63 

Garnett, D. R., Skillman, E. D., Dufour, R. J., Shields, 
G. A., 1997a, ApJ, 481, 174 

Garnett,D.R., Shields,G.A., Peimbert,M., Torres- 
Peimbert, S., Skillman, E.D., Dufour, R. J., Ter- 
levich,E., Terlevich,R.J., 1999, ApJ, 513, 168 

Gonzalez Delgado, R. M., Perez, E., 2000, MNRAS, 317, 
64 

Izotov, Y. I., Thuan, T. X., 1998, ApJ, 500, 188 
Izotov, Y. I., Thuan, T. X., 1999, ApJ, 511, 639 
Izotov, Y. I., Chaffee, F. H., Foltz, C. B., Green, R. F., 

Guseva, N. G., Thuan, T. X., 1999, ApJ, 527, 757 
Kingsburgh, R. L., Barlow, M. J., 1994, MNRAS, 271, 

25 

Liu, X.-W, Storey, P. J., Barlow, M. J., Danziger, I. J., 
Cohen, M., Bryce, M., 2000, MNRAS, 312, 585 

Luridiana, V., Peimbert, M., Leitherer, C, 1999, ApJ, 
572, 110 

Luridiana, V., Peimbert, M., 2001, ApJ, in press 
McKenna,F.C, Keenan,F.P., Kaler,J.B., Wick- 

stead,A.W., Bell,K.L., Aggarwal,K.M., 1996, 

PASP 108, 610 
Mathis, J. S., Torres-Peimbert, S., Peimbert, M., 1998, 

ApJ, 495, 328 

Oey, M. S., Dopita, M., A., Shields, J. C, Smith, R. C, 

2000, ApJS, 128, 511 
Olive, K. A., Steigman, G., Skillman, E. D., 1997, ApJ, 

483, 788 

Osterbrock, D.E., 1989, "Astrophysics of gaseous nebulae 
and active galactic nuclei" 



8 



STASINSKA 



Pagel B. E. J., Edmunds M. G., Blackwell D. E. et al., 

1979, MNRAS, 189, 95 
Peimbert, M., 1967, ApJ, 150, 825 

Pena, M., Stasiriska, G., Esteban, C, Koesterke, L., 

Medina, S., Kingsburgh, R., 1998, A&A, 337, 866 
Pena, M., Stasiriska, G., Medina, S., 2001, A&A in press 
Ratag, M. A., 1991, thesis, Univ. Groningen 
Ratag, M. A., Pottasch, S. R., Dennefeld, M., Menzies, 

J., 1997, A&AS, 126, 297 
Richer, M. G., 1993, ApJ, 415, 240 
Richer, M. G., Mc Call, M. L., 1995, ApJ, 445, 642 
Richer, M. G., Mc Call, M. L., Stasiriska, G., 1998, A&A, 
340, 67 

Rudolph,A.L., Simpson, J. P., Haas,M.R., Erickson,E.F., 

Fich,M., 1997, APJ, 489, 94 
Skillman, E.D., 1989, ApJ, 347, 883 
Simpson,J.P.,Colgan, S. W. J., Rubin, R. H., Erick- 

son,E.F., Haas,M.R., 1995, ApJ, 444, 721 
Spitzer, L., Jr., 1948, ApJ, 107, 6 
Spitzer, L., Jr., 1948, ApJ, 109, 337 
Spitzer, L., Jr., Savedoff, M. P., 1950, ApJ, 111, 593 
Stasiriska, G., 1978, A&AS, 32, 429 
Stasiriska G. 1980, A&A 85, 359 
Stasiriska, G., Leitherer, C, 1996, ApJS, 107, 661 
Stasiriska, C, Richer, M. G., Mc Call, M. L., 1998, A&A, 

336, 667 

Stasiriska, G., Schaerer, D., 1999, A&A, 351, 72 
Stasiriska, G., 2001, in "The evolution of Galaxies. I. 

Observational clues", eds. J.M. Vflchez, G. Stasiriska 

& E. Perez, Kluwer, in press 
Stasiriska, G., Schaerer, D., Leitherer, C, 2001, A&A, in 

press 

Tovmassian, G., Stasiriska, G., Chavushyan, V. H., 
Zharikov, S. V., Gutierrez, C, Prada, F., 2001, A&A, 
in press 



Full addresses go here 



